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ABSTRACT 


This study examines the formation of a front over a 
mountain due to a pre-existing deformation field. The 
dependent variables are assumed to be independent of x. The 
hydrostatic Boussinesq equations are used with no diffusion 
of heat or momentum. A modified coordinate system similar 
to Phillips' sigma system is used. The model is bounded at 
the top by a rigid plane and periodic boundary conditions are 
used in the horizontal. The solutions show the formation of 
a front within a finite period of time that tilts toward the 
cold air. The frontogenesis experiments with a mountain of 


wavelength 3600 km showed negligible topographic effects. 
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I. INTRODUCTION 


Williams and Plotkin (1968) have shown that atmospheric 
frontogenesis occurs for large time limits using eme quasi- 
geostrophic equation. Their results are unrealistic in the 
respect that the frontal zone does not tilt with nm 
comparison of the linear and nonlinear hydrostatic primitive 
equations by Williams (1972) shows that the linear solutions 
are close to the solutions of the quasi-geostrophic equations. 
The nonlinear solutions are more realistic and show the forma- 
tion of frontal discontinuities in a finite period of time. 
Hoskins and Bretherton (1972) showed that dynamic processes 
in addition to horizontal deformation can produce a front. 
Bannon (1983) and (1984) investigated quasi-geostrophic and 
semi-geostrophic frontogenesis over topography.  Monotonically 
sloping topography beneath a deformation leads to the forma- 
tion of a surface discontinuity in the absence of an initial 
horizontal thermal gradient. A front associated with a 
translating deformation field over an isolated orographic 
feature is affected by the topographically forced divergence 
field as the front passes over the moutnain. Convergence 
will increase frontogenesis and divergence will decrease it. 

The purpose of this study is to investigate the effect 
Of topography on a front which arises from deformation. 


Periodic boundary conditions are introduced because they 


TZ 


simplify the mountain flow simulation. The deformation is 
independent of height. The Boussinesq equations are used 

in which the dependent variables are axi. of x. A 
coordinate system similar to that of Phillips (1957) is used. 
The equations are solved numerically and output is shown in the 


form of vertical cross-sections in the y-g plane. 
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II. BASIC EQUATIONS 


The basic equations used in this study are the Boussinesq 
equations which neglect the coi of the atmosphere. 
The hydrostatic approximation is made and heating and friction 
are neglected. These equations will be applied to a model 
in which the flow is bounded above by a rigid horizontal 
plane. 

The Boussinesq equations [Ogura and Phillips, 1962] may 


be expressed in the z-coordinate system as 


> — ^ 
QV > > ov = x E iə 
nc id V:VV + w wm = Vo f kxv , (1) 
9 0 > 00. 
3t + v:VO + w 72 = 0 , (2). 
oz 0 
0 
> OW — 
Vev + 32 = 0 , (4) 
vhere 
cail 
0 E T ) - 0 ¿ (5) 
"Bo 0 
KO 
d = Ge (EPs r (6) 
Ue Po 0 
and 
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In this model, the Coriolis parameter is constant, and the 
only quantity which varies with longitude is 6$. 

To simplify the finite differencing of the equations, a 
modified coordinate system is introduced in which the lower 
boundary surface becomes a coordinate surface. This system 
permits the application of the lower boundary condition at 
the surface of the mountain. The coordinate system is similar 
to the sigma system introduced by Phillips (1957). 

Following Walton (1968) and DeBoer (1970), the vertical 


coordinate, Zeta, is defined as 





Z-Z, 
where D and z, are defined in Fig. l and Z = D-z_. with 
these definitions, ¿ = 0 at the lower boundary and ¢ = 1 at 


Fir upper boundary. 
Dcanscormabren ofteegsations (1) through (4) from the 


z-coordinate system into the z-coordinate system yields 


> > 
AA e SV Z - vo 99/51 me m 
E + v-Vv + E Yala Vo z | E (Vu 2, mBExv., (9) 
90 = O 
E + v:V0 + E зә É Cea (10) 
55... (11) 
dí 9, 
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a 


V. (ZV) + Z a - 0. (12) 
In equations (9) through GER V is the horizontal gradient 
operator on Zeta surfaces. 
It is desirable to put the equations in flux form in order 


to satisfy certain computational properties. Equations (9) 


through (12) become 





9 (M ) 
x d d EC 
E + Sp D + ər (6uz) - Lx T o ; MCN 
9 (M. ) 
ә 9 (Z иб a. 278 
3t + ә + Tala = Z öy + do 250: 15 FM. ; 
(14) 
OH d d - E 
E gy My? * ӱ(092) = О, (15) 
= E 57 28 ! (163 
0 
oM 
where M = uż, M = vZ and H = 02. 
X y 


The boundary conditions are 
s = 0 at c = O , (18) 


and 
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For this study, the flow will be independent of x and periodic 
in y. $ is the only variable allowed to vary in the x- 
direction, so — = 0 for all other variables. 


In order to obtain an expression for $, equation (16) 


may be integrated with respect to z giving 


с 
ф = f -024с̧жс. | (20) 
0 "0 


The constant of integration, C, can be eliminated by taking 
the vertical average of equation (20) and subtracting from 


equation (20) which gives the following 


Fz | 


с 
ф = g-[ ezac-f o zac) + $ . (21) 
0 0 0 


In the above equation,' 'represents the vertical average of 
a given variable. 

A value for $ must be found before a solution for Y can 
be determined. First, take the vertical average of (17) and 


apply the boundary conditions (18) and (19) which gives 


I 
o Гул s 
— M — = 2 — 0 e (227 
oy y = 


This equation shows that M, must be a constant in order to 


satisfy the boundary conditions. Next, take the vertical 
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average of equation (14) and differentiate with respect to 


y which yields 





e 1202 ә ОР. Дәр Se 
mm y ar Tamil = ay Zs + BR (um бә = ЕМ |) . (23) 


The first term in equation (23) is identically zero from (22) 


and the resulting equation becomes 


This equation can be solved diagnostically for ` This equa- 


tion must be applied in order to satisfy the boundary condi- 


mas (18) and (19). 


To obtain an expression for 2, which is a constant, 


equation (21) is differentiated with respect to x which 


ylelds 
E — | E 
dl ə, 90 > 90 90 
x = EU qdo) sx dU * sx - 
0 0 
Since $9 = O in this study, 


and this quantity is assumed to be independent of y. 
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Consider the domain averages of equations (13) and (14) 


which are 


dU E 20 Г=и 
SEU = LC tM (27) 
= a TE a 
d (M) s EMEN ET EE a 
pa = 2 Sy + 9; ^ ki Lzy f Мән ләк (28) 
where — represents the y average. Initially, = is computed 
geostrophically as 
30 Ee 
ox ` z y t=0 a: 


If this relationship is used for all time and if the mass 

flux in the x-direction is zero initially, then M will be 

zero for all time. The first two terms on the right side 

of equation (28) represent the mountain drag. The sum of 

these two terms for mid-latitude flow has been shown to be 
negative in the mean. This will lead to a continual loss of 
y-momentum. This kind of loss would make a quasi-steady final 
state impossible. An artificial source of momentum is required 


to maintain the mean meridional flow. This artificial source 


ra 


of momentum is provided by holding — constant as obtained 


in equation (29). In this case, if M, is decreased below 


rd i Pe ui : 
the initial geostrophic value, M. will increase and restore 


the momentum loss. Thus, in the long term mean, the mountain 
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terms are balanced by the mean Coriolis force. However, for 
the horizontal scales to be considered in this study, this 
momentum source may not be needed since — scale moun- 
tains cause very little wave drag. 


The dependent variables are to be divided as follows: 


u(x,y,o,t) - Uq(x,y) * u(y,t,t) ; 

yOu) = Va Gc y) uk”), 

L(y,t,t) = çt(y,c,t) , (30) 
CRG ct) IAN r И. 

ME) =D IO AE?) 


The basic wind deformation field is chosen to be periodic in 
y and to have zero vorticity and divergence. Thus, the 
deformation field is given by 


О (х,у) = –О Lsinh (7) sin (2) , (31) 


d 


b" 


Vą (XY) -D4 Lcosh (7) cos (+) , (50) 


EF 


vhere Da and L are constants. This deformation field was 


used by Stone (1966) in his study of quasi-geostrophic 


SCH 


frontogenesis. Since the flow is assumed to be independent 


of x in this study, equations (Sl) and mamme evaluated at 


x = 0 resulting in 
U,(0,y) = 0, ( 3:99 
- - Y 
V4 (Oy) Da Pcosİir) . (34) 


Substituting the expressions in equations (30) through (34) 
into equations (13) through (15) and droppine the prime 


notation yields 


o (V-Z) 


du du d d o El 
2 3t + .. “öy "U DE + = a z (UZ) 
= - ZAŚ + vz, —. 
= 
o (V-Z) | 
öv 9 d öv d 
2 3t + mine + v 5 m Va ә az (bua) 
- 0557707 19 94 . 36 
Z öy + бу 982(2-1) m ЈИ, (36) 
90 9 d P 
9v oz 2. = 37 
Z 5 + . + z (507) 0, (3% 
where Va = -D; Lcos (7). This last set of equations is the 
set that needs to be solved for u, v and 8. In this study, 


equation (24) does not have to be solved if M, is restored 


to the initial constant value at every timestep. 
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III. FINITE DIFFERENCE FORMULATIONS 


Staggered grids are used in both the horizontal and 
Rertical. The vertical grid is divided into K-l levels of 
constant Az (Fig. 2). Quantities defined for the layers will 
be denoted by integer subscripts and the z coordinate value 


will be denoted by zl The quantities 0, u, v and ó are all 


e 
defined on the integer levels. The vertical motion, 2, on 
the other hand, is defined on the top, bottom and interface 
levels. These levels will be denoted with half-integer 


Subscripts and the vertical coordinate will be identified by 


C41 /2° The horizontal grid is divided into J-l sections of 


constant resolution (Fig. 3). Momentum variables are defined 
on the integer subscripted points denoted by E The momentum 
variables include u, v, Mer RW and Va The mass variables, 


WO, Z and + are defined on the half-integer subscripted 


points. These points will be denoted as YI Values 


uz”: 


needed off their respective grid points vill be identified 


by a caret and vill be defined as 


mür. 
Y 


ro 
14.172] 


(Y. 11721) 8541172 
J ии 


U > 
Il 


, (38) 


for quantities defined on the half-integer points and as 
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ÓK1+1/2 d 
Am cuve 


Gk1-1/2 9 G 


ÓWETZ Ne 
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Ç 3⁄2 d 
Gale n uU, V, d p 


61/2 әс̧ 


Fig. 2. ТҺе vertical grid, showing the indexing 


convention and staggering of the variables 
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E _ Pay p n  — j /5)S. A 
ua 217. raq à - 
Jj+1 j 
for quantities at the integer points. Periodic boundary 


conditions are used in this study for adii quantities. 

Equations (11), (13) through (15) and (17) are solved by 
using finite differences in y, 7 and t. The finite differen- 
cing scheme is chosen to conserve the sum of kinetic energy 
and potential energy. 

The time differencing scheme uses a leapfrog scheme with 
a periodic use of the Euler backward scheme. The Euler 
backward scheme is used to eliminate the computational mode 
arising with the leapfrog scheme [Haltiner and Williams, 
1980). 

Equations (11), (17) and (13) through (15) may be written 


in finite difference form as 


M 
u.c bu 


M M ` M 
Y "7 u. UL y=", SAMA ._ 


JOSE Tur - E2 
(22) Job "Au. ve (LK Me kel, 
2 JASTA 2 j,k-1/2 2 
>kH/2 TAS 
= Em 5 2 (40) 


26 


M M M 
öv Qna Ok HIEL ab 9 (Ypk Yjlgk "ik jk 
әв —— —— —— S———— ——  — 
: A 1-172 


y NS +V. be 
: De ae jk Vj,k-1 


“ET” ^k-1/2 


EN 


5141/2 23-172 


16 UZ 1/2 


. ф. m. RA 
=- -7 J+1/2,k 'j-1/2,k iL 2-2), (2-2) 


J Hay "0 


^^ 


" x k 23 k ' (41) 


s M (j+3/2,k j+1/2,k, Мм Qjd/2 1/2, 
: ДМК 2101,К 2 Уј,к 2 


RE 2 dE EM) = Y. 


(22) иә јл/2,Ҝ) (с2) C372 x 372,1, 
4-7 j+1/2,k+1/2 2 Eer > 


ökal/2 7 ök-1/2 


- Ü, (42) 


^jil/2,k*l ^ ?j«1/2,k Gs gur I 568424 k 
ит. 2 


ECL Cy 0 
(43) 


and 


M =M : ‘ 
y a 

j+l,k *j,k _ -7 A a2 kE 
YB AWZ | ><) 


5 +1 m mm 7 Ök-1/2 
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IVe wEE5UPES 


All numerical results to be shown use the following 


numerical values for the constants: 


a = 12 K 
B. 375 m 
D = 9 km 
Da - Ro sec + 
f a (45) 
g = 9.81 m/sec 
L = 3600 km 
99 = 300 K 
"Us - 4 K/km 
öz 


The folloving finite-difference increments are used. 


At = 180 sec 
53257 (46) 
Ди – 122, 


A. MOUNTAIN WITH NO DEFORMATION FLOW 

The initial conditions are obtained from the quasi- 
geostrophic equations for steady state flow over topography 
[Walton, 1968; DeBoer, 1970]. The initial potential 


temperature field is given by 
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90 
R s sinh[a(z-D) J p aa y 

AR 2,0) = J + B EI кист COS (7) = cos (7) : 77. 

99. 
The quantity E is the constant initial lapse rate and 

1757 о 
. Qu. xs 
m... FL (8 52 : on 


0 


Substituting equation (47) into the thermal wind equation 


yields 
ou . ig lude sinh [a(z-D)] 515 (Y) - 25 sin (291 (49) 
SE t 09 Ы öz sinh (aD) L T L i 


This equation can be integrated with respect to z to give 


90. 1/2 

- 6 coshlo (z-D)1 ,Y 

u(y,z,0) MER T sinh (aD) sin (7) 
- 22 sin (2) (Z -3) (50) 


The initial v field can be shown to be 


cosh [a ( z-D) ] 


MZ) = V + GVB SEENEN 


y 
cos (r) d (51) 
The initial w and à are obtained from equations (4) and (3). 
It is necessary to rewrite the initial conditions in the 
zeta-coordinate system. From equation (6), equations (47), 


(50), and (51) can be rewritten as 


— 


985 sinh[aZ(t-1)] 


= 1 ə _2a 1 na 
on S a 0) c B TE Sim sin (7) Ñ sin (7) , GAR 
a qua a 
i g S coshilo ZA. y 
u(y,c,0) ar SE sinh aD) nn 
-2a sin (EZ) (¢ -5)D (53) 
T m e i 
and 
» cosa sul y 
121700) V + a VB EXE cos (+) š (54) 
These initial conditions are shown in Figs. 4-7. 1 and $ are 


given by equations (10) and (9). 


The value of 0, Was computed using the formula 
98 
2 Sri - 
Es = ol? dy. (339 


The profile of the mountain used in this study is given by 


| (56) 


вн 


2. = E cos 


and is shown in Fig. 8. 

The initial conditions show anticyelonic voxrticiCYWGNNJE 
the ridge and cyclonic vorticity over the valley. There is 
divergence on the upslope and convergence on the dovnslope. 


These perturbations die off vith height above the surface. 
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(W) 


Figs. 9-12 show the results at t = /2 houes.  SITNEGNENE 
original fields were analytic, only a slight adjustment to 


the mountain can be seen in the v, and hence, .the 2 fields. 


B. DEFORMATION OVER FLAT TOPOGRAPHY 


For this case, the initial potential temperature field 


contains a horizontal variation. The initial field is given 
by 
90 
- S (z -D) _ 2a y 
08(y,z,0) = 37 (2 5 - cos (r) : (52) 


The initial u field is again derived using the thermal vind 
relationship and can be shovn to be 
= = ии D 
im — sin (F) (z -5) . (58) 
The initial v field is identically zero and thus, w is fODGNE 
to be zero by equation (3). 
It is again necessary to write equations (48) and (49) 


in zeta-coordinates as 


90 


9(у,с,0) = —ӱ8(с-1)0 - “ZP cos (D) , (59) 
and 
u(y,0,0) = - ŚBsin($) (z -5)D . (60) 
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v and с are zero initially and qis-given by equation (9). 
amas. 13 and 14 show the initial conditions and Fig. 15. 
shows the deformation flow given by equation (34). 

The front is developed by 24 hours as shown IES. 16-19. 
The frontal discontinuity forms along the point where v = 0 
and tilts toward the colder air. Cyclonic shear can be seen 
along the frontal boundary. The relative vorticity is less 
than zero on the north side of the front and is greater than 
zero to the south. 

By 48 hours, Figs. 20-23 show the strengthening of the 
frontal discontinuity. The tilt remains the same as at the 
24 hour point. On the other side of the domain, the tempera- 
ture gradient is weakening, which indicates frontolysis 


S@eurring in that region. 


C. MOUNTAIN WITH NON-MOVING DEFORMATION 
The previous two cases are combined in order to investi- 
gate the combined effects of a deformation flow and a 


mountain. The initial fields are given by 


6 . 
B s nl ace m y 2а y 
2157 erg, DI ә Ee T B 37 sinh (aD) sı (7) TT COS (F) , (61) 
Ig 12 
b a S cosh [az (c-1)] хи 
Uyam» = dn 37) AA sin (r) 
ú 2a . y 
— Sin (7) г (62) 
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ha cosh[aZ2(zg-1)] y 
V(y,5,0) = V + aVB EOXENEGD) Cos (p) 7 5) 
and š is given by equation (10). “These fields can be seen 


in Figs. 24, 25, 26 and 27. The deformation is the same as 
in the previous case. 

As in the previous case, a tightening of the temperature 
gradient has occurred by the 24 hour point, although not as 
strong (Figs. 26-29). The front has moved farther to the 
north and the tilt is smaller than in Case B. Throughout 
the next 24 hours, the front continues to move to the north 
and passes over the mountain ridge 157700. The 
gradient veakens as the front passes over the ridge and the 
tilt is negligible. In this experiment, the front weakens 
because it is advected away from the frontogenetic part of 


the deformation field. 


D. MOUNTAIN WITH MOVING DEFORMATION 

The next four cases to be considered start with the same 
initial conditions as in Case C (Figs. 24, 25, 26 and 27). 
However, in these cases, the deformation moves to the north 
at a speed of 10 m/sec. Each case examines a different 
starting place for the maximum northward deformation flow. 
The cases have the maximum flow over the valley, on the 
upslope side of the mountain, directly over the ridge, and 


on the downslope side of the mountain. 
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1. Maximum Deformation over the valley 


The initial deformation is shown in Fig. 15.: NUM 
the previous cases, the front has formed by 24 hours (Figs. 
34-37). There is a small tilt toward the colder air in the 
north and the front has moved farther up the mountain than 
in Case C. The upper level front has passed over the moun- 
tainfğ r idae: 

As the deformation flow continues to move northward, 
so does the front. During the next 24 hours, the speed of 
the front slows down, but the gradient in the temperature 
field tightens (Figs. 38-41) and forms a discontinuity. The 
cyclonic shear around the front has increased and there is 
strong vertical motion in that region. 

A careful comparison of the potential temperatures 
at t = 24 hours between Figs. 24 and 34 shows very little 
change between the flat experiment and the mountain experi- 
ment, except for the translation in Fig. 34. The same con- 
clusion holds at 48 hours. Bannon (1983,1984) has shown that 
mountain effects on frontogenesis will be important when the 
mountain forced convergence is of the same order as the 
frontal convergence. Comparison of v in Figs. 6 and 18 shows 
that the mountain forced convergence is only about 25 percent 
of the mountain forced convergence at 24 hours. This ratio 
is even smaller later in the frontal evolution. With the 
large scale mountain used in this and later experiments, the 


effects of the topography are unimportant. 
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2. EE sShifted One-Quarter Wavelength 


In this experiment, the deformation field. is shifted 
one-quarter of a wavelencth to the north with respect to 
the initial fields (Fig. 42). The topography is also unim- 
portant in this experiment. At t = 24 hours, there is only 
a slight indication of frontal formation in the temperature 
field, but a much stronger suggestion of the front in the u 
field. The u field, however, shows little tilt (Figs. 43- 
46). Williams (1968) has shown that this phasing of the 
deformation field is inefficient for frontogenesis. 

3.  Reversed Deformation Field 

As is shown in Fig. 47, the initial deformation flow 
is in opposite phase with the original deformation field. 
The front can be seen at 24 hours in Figs. 48-51. As with 
Previous cases, the front tilts toward the cold air, but in 
this case, the tilt is toward the south. This is a result 
of the reversed deformation field. The fields are mirror 
images Of the fields in the basic experiment except for the 
Shift. 

4.  Deformation Reverse of Case 2 

The last case to be considered is the one in which 
the maximum deformation flow is the reverse of Case 2 (Fig. 
52). The first indication of frontal activity can be seen 
24 hours later (Figs. 53-56) especially in the u field. 
Strong cyclonic activity can be seen here with no evidence 


of a tilt. The temperature field shows only a small tightening 
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cs r from the initial conditions. This is very 
similar to Case 2, but the maximum shear is at the lover 


boundary. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


In this study, a numerical model of frontogenesis was 
constructed which included — topography. Periodic 
boundary conditions were included and a periodic deformation 
field was introduced. The model was successfully tested with 
frontogenesis with no topography and for stratified flow over 
a mountain with no deformation field. Combined experiments 
with topography and a moving deformation field were carried 
out. The experiments showed negligible effects on fronto- 
genesis by the topography. 

The next step in this research is to introduce smaller 
scale topography which should lead to a noticeable effect 
as suggested by the frontogenesis theories of Bannon (1983, 
1984). Other effects including vertical wind shear could 


be included later. 
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